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We introduce a partial state fidelity approach to quantum phase transitions. We consider a 
superconducting lattice with a magnetic impurity inserted at its centre, and look at the fidelity 
between partial (either one-site or two-site) quantum states. In the vicinity of the point of the 
6J[)| quantum phase transition, we observe a sudden drop of the fidelity between two one-site partial 

states corresponding to the impurity location and its close vicinity. In the case of two-site states, 
• the fidelity reveals the transition point as long as one of the two electron sites is located at the 

impurity, while the other lies elsewhere in the lattice. We also determine the Uhlmann mixed 
state geometric phase, recently introduced in the study of the structural change of the system state 
eigenvectors in the vicinity of the lines of thermal phase transitions, and find it to be trivial, both 
for one- and two-site partial states, except when an electron site is at the impurity. This means 
that the system partial state eigenvectors do not contribute significantly to the enhanced state 
distinguishability around the point of this quantum phase transition. Finally, we use the fidelity to 
analyze the total amount of correlations contained within a composite system, showing that, even 
for the smallest two-site states, it features an abrupt quantitative change in the vicinity of the point 
. of the quantum phase transition. 

PACS numbers: 03.67.-a 05.70.Fh 75.40.Cx 



Recently, considerable developments have been made within the research born out of the collaboration between 
physicists from two rather distinct fields, viz those of quantum information and computation and those of condensed 
matter and many-body physics. Most of this interdisciplinary work has focused on the realization of quantum 
information protocols with real many-body systems (for an overview, see [lj|). Yet, a considerable development 
• has been accomplished by moving in the "opposite" direction, ie applying concepts and techniques widely used within 
the field of quantum information and computation to the study of problems relevant within the field of condensed 
matter physics. One of such topics, recently researched to a great extent, is that of the use of entanglement in the 
\ study of zero-temperature quantum phase transitions Q , as well as general thermal phase transitions (for an overview, 
■ see [3(, chapters 19-21, and references therein). 

One application of the quantum information techniques to problems dealing with macroscopic many-body systems 
is the use of the (ground) state fidelity Q in studying (quantum) phase transitions. It was first reported in [f|, 
where the sudden drop of the fidelity along the regions of quantum phase transitions was observed in the Dicke 
and the XY models. Consequently, this method was applied to the cases of free Fermi systems Bose-Hubbard 
models 0, H[, spin chains [ij [l(| and in connecting the fidelity with the renormalization group flows 11]. Also, 



the cases of non-symmetry breaking quantum phase transitions that are not of the Landau-Ginzburg- Wilson type 



were studied [12], Il3l . Il4j . A more formal and general study based on differential geometric insights, establishing a 



connection between quantum phase transitions and the relevant thermodynamic response functions, was developed 



in [15|, [16[ (see also [iZ[). Finally, using the mixed state fidelity, and extending the parameter space to include the 



temperature, quantum phase transitions were studied in 18], while the thermal phase transitions were explored in 
[TBL Il7l. Il9j|. It is also worth noting here that the recently studied connection between the quantum phase transitions 



and Berry geometric phases [2Q| was established in [15|, |16| on a more formal level, unifying the fidelity and the Berry 
phase approaches to quantum phase transitions. Subsequently, the Uhlmann mixed state geometric phase [2ll ] (the 
generalization of the pure state Berry phase) was introduced in [r| as a tool in exploring the change of the system 
thermal equilibrium state eigenvectors in the vicinity of the line of a phase transition. 

Given two physical systems a and b, their quantum states, being either pure or mixed, are given by density operators 
p a and pb- Since according to quantum mechanics it is in general not possible to unambiguously distinguish with 
probability one between two given states [H, various measures of quantum state distinguishability have been 
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introduced (for an overview, see |22j). One of the most widely used measures is (quantum) fidelity [4|, given by the 
expression 



F{Pa,Pb) = Tr 



' PaPby 



(1) 



The numerical value of the fidelity ranges from F = 0, in the case of fully distinguishable states [37(, to F = 1, when 
the two quantum states are completely indistinguishable from each other (ie they are identical) . As it happened with 
most other quantum state distinguishability measures, the introduction of the quantum fidelity F was also motivated 
by the analogous classical quantity. Classical fidelity F c , a measure of distinguishability between two probability 
distributions {p a {i)} and {pb(i)}, is defined as F c (p a ,pt,) — J2i y/Pa(i)Pb{i) ■ Measuring an observable A in states 
p a and pb, one obtains probability distributions {p^ (i)} and {p^(«)}) respectively, whose mutual distinguishability 
is given by the classical fidelity F c (p^,p^). It can be proven (see [22j]) that for any observable A, the inequality 
F(p a ,Pb) < F^p^Pb) holds, and there always exists an optimal observable A op for which the equality is achieved. 
In other words, quantum fidelity F(p a ,Pb) gives the optimal value for the distinguishability between two quantum 
states, when compared through probability distributions {p^ («)}) and {p^(i)} (note though that it does not give an 



optimal observable A op ). Using the modulus of an operator R, \R\ = (RR 1 ) 
for the fidelity fl: 



1/2 



one obtains an alternative expression 



F(p a ,Pb) = Tr 



= Tr 



)t = Tr 



' PaPb^ 



(2) 



The fidelity approach to quantum phase transitions is based on the use of the above quantity ([I]) in distinguishing 
macroscopic quantum states of many-body quantum systems. Consider a system given by the set of Hamiltonians 
H(q), where q is a parameter representing a set of interaction coupling constants, which at T = is in a pure ground 
state |$) = \®{q)). Then, the sudden drop of the fidelity Q |$(g))($(g)|) = between 

two ground states |$(<?)) and \<&{q)) obtained for two close parameter points q and q = q + 5q, with 5q being small, 
can indicate the points of quantum phase transitions. The basic logic behind this idea is strikingly simple — two 
quantum states defining different macroscopic phases are expected to have enhanced distinguishability that would 
quantitatively exceed that taken between the states from the same phase. In our study, instead of the global ground 
states that describe the overall system, we will study partial quantum states of certain subsystems of the system. If the 
overall system S is divided into two parts A and B, then we will study the mixed state fidelity F(pA(q), Pa{q))i where 
PA (q) = Tr# | $(<?)) ($>(q) |, and analogously for Pa(q)- We will show that in the system considered, a conventional BCS 
superconductor with an inserted magnetic impurity, a sudden drop of the above mixed state fidelity can define the 
point at which the (first order) quantum phase transition occurs, even when applied to the smallest subsystems of 
the overall system. 

THE MODEL 

In this work we consider the quantum phase transition induced by a magnetic impurity inserted in a conventional 
superconductor 23, [iT 



25, 26, 27, 28, 2 



Q|. Consider a classical spin immersed in a two-dimensional s-wave 
conventional superconductor. We use a lattice description of the system. In the centre of the system, i = l c = {x Cl y c ), 
we place a classical spin parametrized by 



S 



cos ipe x + sin ipe z . 



(3) 



The Hamiltonian of the system is given by: 

(ij)(T ity i utj' 

where the first term describes the hopping of electrons between different sites on the lattice, ef is the chemical potential 
(equal to the Fermi energy at T = 0), the third term is the superconducting s-pairing with the site-dependent order 
parameter Aj, and the last term with J > is the exchange interaction of an electron at site i = l c with the magnetic 
impurity. The hopping matrix is given by = tSj^+s where S is a vector to a nearest-neighbor site. Note that both 
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the indices I and i,j G {1, 2, . . . N} specify sites on a two-dimensional system (N is the number of sites). We will take 
energy units in terms of t (t = 1), and ep = — 1- 

Since the impurity spin acts like a local magnetic field the electronic spin density will align along the local spin. 
For small values of the coupling there is a negative spin density around the impurity site. At the impurity site it is 
positive, as expected. For larger couplings the spin density in the vicinity of the impurity site is positive. At small 
couplings the many-body system screens the effect induced by the impurity inducing fluctuations that compensate the 
effect of the local field in a way that the overall magnetization vanishes. However, for a strong enough coupling the 
many-body system becomes magnetized in a discontinuous fashion. One interpretation is that if J is strong enough 
the impurity breaks a Cooper pair and captures one of the electrons, leaving the other electron unpaired, and thus the 
overall electronic system becomes polarized. The impurity induces a pair of bound states inside the superconducting 
energy gap, one at positive energy (with respect to the chemical potential), and another at a symmetric negative 
energy. Even though the spectrum is symmetric the spectral weights of the two energy levels are not the same and 
their spin content is also distinct. The analysis of the local density of states (LDOS) 27, HI shows that for small 
coupling the lowest positive energy level has only a contribution from spin \ and the first level with negative energy 
(symmetric to the other level) has only a contribution from spin J.. The magnitude of the spectral weight at the 
impurity site is different for the two states. Considering a higher value for the coupling one finds that the levels 
inside the gap approach the Fermi level. There is a critical value of the coupling for which the two levels cross in a 
discontinuous way such that it coincides with the emergence of a finite overall magnetization. After the level crossing 
occurs, the nature of the states changes. The positive energy bound state has now only a contribution from the spin 
I component and vice-versa, the first negative energy state has only contribution from the spin component | — hence 
as the level crossing takes place the spin content changes. 

The diagonalization of this Hamiltonian is performed using the Bogoliubov-Valatin transformation in the form: 



C»T = ^[Unii^^n- V*(i, T)7n] 
n 



(5) 



Here n is a complete set of states, u n and v n are related to the eigenfunctions of the Hamiltonian Q, and the new 
fermionic operators j n are the quasiparticle operators. These are chosen such that in terms of the new operators: 



Trillin, 



(6) 



where E g is the ground state energy and e n are the excitation energies. As a consequence: 



H ^ 

A± i In 



t 



(7) 



The coefficients u n (i, a), v n (i, a) can be obtained solving the Bogoliubov - de Gennes (BdG) equations [32| . Defining 
the vector: 



( Un(i,1) \ 
Vn(i,l) 

V W "(*>T) / 



(8) 



the BdG equations can be written as: 



Htp n (i) = e n ip n {i), 



(9) 



where the TL matrix at site i is given by: 



H 



/ —h — Ep — J sin <p 
A* 

— J COS ifi 

V o 



A, — J cos tp 

h + ep — J sin ip 

— h — ep + Jsmip 

— Jcosip A* 





— J COS (fi 

A, 

h + ep + J sin tp J 



(10) 
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where h — ts$ with sgf(i) = f(i + 5). The solution of these equations gives both energy eigenvalues and eigenstates. 
The problem involves the diagonalization of a (47V) x (47V) matrix. The solution of the BdG equations is performed 
self-consistently imposing at each iteration that 

Ai = ^[(cirCil) - (cijCi T )], (11) 

where g is the effective attractive interaction between the electrons. Using the canonical transformation this can be 
written as 

A» = -g ^2 |/n(e«)K(«,T)<(i,l) +Un(i, l)v*(i, t)] - ~[u n (i, ?)<(«, I) +u n (i, l)v*(i, t)]| (12) 
where lod is the Debye frequency, and f n (e n ) is the Fermi function defined as usual like 

= (13) 

where T is the temperature. We note that the Bogoliubov - de Gennes equations are invariant under the substitutions 

e n -» -e n , u(T) - w(T), w(T) - «(T), «U) - «U) - -«U)- 



RESULTS 



As we have already noted, in this paper we focus on the fidelity between the partial states rather than the overall 
ground states, as it was the case in the previous studies 

system is in the pure ground state |$) and has a density operator p = |$)($|. Thus, if we divide the whole system 
in two subsystems, say A and B, then the partial mixed state, given by the reduced density operator pa for the 
subsystem A, is defined as 



p A = Tr B p, 



(14) 



where Ttb[-] represents the partial trace evaluated over the Hilbert space TLb of the subsystem B. In particular, we 
will calculate the matrix elements of one-site and two-site partial mixed states, given in general for different sites i and 
respectively, and different values of the coupling constant J. In order to infer the point of the quantum phase 
transition, we will calculate the fidelity between two one-site and two two-site states, given for two close parameters 
J and J + S J. Further, to better understand the influence of the quantum phase transition point on the structure 
of partial states (and therefore, to a certain extent, the overall state as well), we examine the fidelity between two 
different partial states, given for the same values of the coupling constant J, but different sites, for both cases of one- 
and two-site states. 

In many-body systems a description in terms of quantum pure states and wave functions as their representatives 
in a certain basis (e.g. position, momentum, energy,...) is quite involved and the second quantization is the natural 
way to perform any calculation. The mixed states are however easily defined in terms of pure states and their matrix 
elements in some basis, and this approach is indeed not difficult to implement using Fock states. The matrix elements 
of the density matrix are simply defined in terms of correlation functions of the whole system. For instance, in the 
case of the single-site partial mixed states it can be shown that, using local basis states B — {|0), 1 1|) | T)> 1 1)}: which 
denote the four possible states — unoccupied, double occupied, single occupied with an electron with spin up and 
single occupied with an electron with spin down, respectively — the corresponding density matrix reads as 



V 



n T )(l 
(c;c T ) 







(44) 






(n T (l-nx)) (ctc T ) 



<4 6 j.) 



(15) 



((l-n T )n|) J 



where the index i denotes the site. The spin and the charge parts decouple. The spin part couples the two spin 
orientations (single occupied states) and the charge part couples the empty and doubly occupied states. The diagonal 
terms of the matrix describe the number of empty sites, the number of doubly occupied sites, the number of spin up 
sites and the number of spin down sites, respectively. The sum of the diagonal terms is equal to 1 due to normalization. 
The matrix is easily diagonalized and the fidelity (|TJ) between two different one-site states obtained straightforwardly. 
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FIG. 1: One-site fidelity between two states taken for two close parameter points J and J + SJ (SJ = 0.05), for the site at the 
impurity (black line), a nearest neighbor (green line) and a site located in the bulk (red line). Note the sudden drop of the 
fidelity for the impurity around the point of the quantum phase transition. 



The correlation functions are easily solved using the representation of the electron operators in terms of the BdG 
quasiparticle operators. Specifically, we can write the single-site density matrix as 

P = ^2pnm\n)(m\ (16) 

n, m 

where \n), \m) £ B are the four basis states described above. Consider now an operator O defined in terms of creation 
and annihilation operators of the electrons, and evaluate 

Tv{6p} =J2(n\d\m)Pmn. (17) 

n,m 

To determine the matrix element p nm of the density matrix it is enough, by inspection, to find which operator O 
is such that only the matrix element (n\0\m) is nonzero. This will tell us which correlation function gives which 
matrix element of the density matrix. Alternatively, one can express the operator |m)(n| in terms of the creation and 
annihilation operators. The results for pa are given by equation (j 1 5[) above. 

These matrix elements are defined in the subspace of the states of one site. However, since we have integrated out 
all the other sites, we may now replace the matrix elements by the matrix elements over the entire system (direct 
product with the states of the remaining N — 1 sites) . 

Computing the fidelity between the two-site partial states is more involved. Given two sites, i and j, the two-site 
basis states are given by the direct product of the one-site basis sets that were defined earlier, and the reduced density 
matrix is now a square matrix of order 16. However, exactly half of the 256 matrix elements are null — the nonzero 
ones being those that correspond to either even-even or odd-odd states, where the parity reflects the number of 
electrons in the sites i and j, respectively. Hence those two sets of states decouple and the reduced density matrix 
can be written as a block diagonal matrix where each submatrix is a square matrix or order 8. Still, the calculation 
is rather complex and for that reason we used the computer program FORM [33| to obtain symbolic expressions for 
the matrix elements involved. Note that one of the correlation functions involves the double occupancies at the two 
sites, and thus one should deal with products of up to eight electron operators (where each electron operator is given 
by a sum of two quasiparticle operators); also, the signs that follow from the anti-commutation relations of fermionic 
operators must be taken into account. 
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FIG. 2: The charge and spin contributions of the one-site fidelity as a function of the parameter J, with 8 J = 0.05, for the site 
at the impurity (upper figure) and a site located in the bulk (lower figure). Note the discontinuity for both spin and charge 
contributions around the point of the quantum phase transition, for the site located at the impurity. 



Fidelity Between Single-site States 

In the following, we present the numerical results for the fidelity between one-site partial states, for a two dimensional 
lattice of size 15 x 15 sites. For computational simplicity, we perform our calculations for this system size. Since 
the effect of the impurity is short-ranged, changing the size of the system only affects slightly the results. First, we 
evaluate the fidelity (fTJ) between two one-site states given for two sites i and j and two close values J and J + S J of 
the exchange parameter, with 5 J = 0.05, such that p a = p(J;i) and pb = p{J + SJ;j). 

For i = j = l c , when the electrons are located at the site of the magnetic impurity, the results for the one-site 
fidelity F\{p a , fib), as a function of J, are presented in Fig. [1] We observe a clear drop, from the otherwise common 
value Fi(p a , pb) ~ 1, for the critical value Jo ~ 1.9 [3l| of the parameter J. As we move one of the electrons from the 
impurity, the drop of the fidelity F\ around the point of the quantum phase transition becomes much smaller: for the 
first neighbors of the impurity, it is still visible, while already for the second-neighbors and the far away bulk sites, 
it becomes negligible. We see that even for the smallest subsystems — one-site spins in the lattice — the fidelity 
can be dramatically influenced by the structure of their states, as long as they are close to the impurity. This is 
intuitively easy to accept since it was shown in the previous work (3ll . [34| , through investigation of various relevant 
physical quantities, such as the gap, the electron density or the total and quantum correlations (entanglement), that 
the impurity affects one-site properties of the system mainly in its vicinity (although it was also shown that for some 
two-site entanglement measures the effects persisted for larger distances). 

A more detailed picture of the fidelity between the one-site states is presented in Fig. [21 where the two distinct 
contributions to the overall change of partial states is presented. The charge and the spin parts of the one-particle 
state (fT5)) lead to two different contributions to the state distinguishability, and we see that they both exhibit a 
discontinuity in the vicinity of the point of the quantum phase transition. 



Further, we have evaluated the Uhlmann mixed state geometric phase 21|, a generalization of the Berry geomet- 
ric phase to mixed quantum states. While the connection between Berry phases, pure state fidelity and quantum 
phase transitions was studied before [TBI . [T5 , [2fj| . the Uhlmann mixed state geometric phase was recently intro- 
duced [lj| to analyze the structural change of the system eigenvectors, with respect to the change of the relevant 
parameter driving the system Hamiltonian (ie J in our case) . The Uhlmann geometric phase 4>g (and the correspond- 
ing holonomy) is determined by the unitary operator V a b, given by the polar decomposition (see for example [l|) 
^Tpa^fpb = \ Vp~aVpb \V a b, where |-R| = (RR)) 1 / 2 represents the modulus of the operator R (for details, see [IBj]). Since 
the fidelity can be written as F(p a ,pb) = Tr\^/p\^/pb\, see equation ([2]), a nonzero difference H(p a ,pb) — F(p a ,pb), 
where H(p a , Pb) — ^[V^aVpb], would imply V a b ^ lab and therefore the emergence of a non-trivial Uhlmann geomet- 
ric phase (j} g ^ 0. On the other hand, for the case of mutually commuting Hamiltonians the system eigenvectors do 
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not change with the driving parameter and we have F(p a ,pb) = Trfy^paVPb ]■ In other words, F(p a ,pb) — H(p a ,pb) 
and the Uhlmann geometric phase becomes trivial, ie 4> g = 0. 

For the one-site states, we find that the Uhlmann geometric phase is always trivial for every site neighboring 
the impurity or in the bulk: the one-site eigenvectors do not contribute significantly to the change of the overall 
one-site partial states at the point of the quantum phase transition. Still, for the impurity site itself, around the 
point of the phase transition Jo we observe a small deviation from zero for the difference H± — F\, which indicates 
a slight "rotation" of the electron's eigenbasis. This result can be qualitatively explained noticing that in the model 
we are considering the third term in the Hamiltonian (j4|), giving the superconducting features, and the last term 
in the Hamiltonian, giving the coupling between the impurity and the lattice electrons, commute with each other 
[38j |. Moreover, for a fixed value of the gap A», changing the impurity-lattice coupling constant J, without changing 
the impurity spin orientation angle ip, will result in no change of the overall Hamiltonian's eigenbasis — only the 
eigenvalues are affected in this case and the system exhibits a first-order quantum phase transition as a consequence of 
the explicit (not avoided) level crossing. Thus, one might expect that also the partial subsystem states would feature 
little change of local eigenbasis, even in the vicinity of the point of the quantum phase transition. This is precisely 
the case in the model considered, apart for the gap A; c given for the electron located at the impurity (see [3l||): for 
every site i, the gap Aj is a slowly changing function of J, while for the impurity A; c exhibits a tt shift and becomes 
negative. Such a behavior of A/ c does affect the difference Hi — F\, but the effect is an order of magnitude smaller 
than for the fidelity. Further, it affects the impurity site only and thus does not lead to the superconductor-normal 
metal type of phase transition. 

Finally, to further investigate the structure of one-site partial states, we have evaluated the fidelity between the 
states for two different sites, and the same value of the parameter J. The results are presented in Fig. [3l Again, 
we see that apart from the state of the site l c located at the impurity, and those at its vicinity (first neighbors), all 
the other one-site partial bulk states are almost the same, for a given common value of the parameter J. In Fig. 
[3] we have presented the results for two distinct values of the parameter, J~ = 1.5 < Jo and J + = 2.5 > Jo, but 
the qualitative picture stays the same for every J~ < Jo and J + > Jq\ the numerical value for the fidelity between 
the bulk state (say, p±) and the impurity state (say, pi) is significantly greater for J~ than for J + , in other words, 
F^(p a ,pb) > F^ (p a , Pb), thus revealing the point of the quantum phase transition. 

Note that when i = l c and J < Jq, besides the peak at the impurity site, in its vicinity the fidelity F\ is smaller than 
the background value, while for J > Jo it is larger. This is directly connected to the existence of oscillations of the 
magnetization around the impurity site (see [3lj| ) . Within the bounds set by quantum mechanics, two quantum states 
can be distinguished by measuring different observables and comparing the outcomes obtained. Magnetization is just 
one of them, and in this case we see that it is optimal, or at least nearly so: for J < Jo, it is easier to distinguish the 
impurity site from a site in its immediate vicinity than from other, more distant sites, simply because the difference 
in the expected magnetization between two sites is bigger in the case of impurity-neighborhood than in other cases. 
Therefore, in this case the fidelity has a drop in the neighborhood of the impurity, ie the distinguishability is larger 
there than elsewhere. A similar analysis holds for J > Jq. 



Fidelity Between Two-site States 

We have also studied the behavior of the fidelity taken between two-site partial states. We considered two distinct 
cases. First, we calculated the fidelity between two-site states for which one of the sites i,j is at the impurity, say 
i = l c , and the other is in its vicinity, j '• = l c + 51 (51 = 1,2,.. .6), for two close parameter values J and J + 5 J, 
with 5 J = 0.05. In Fig. [4] we present the results for the two-site fidelity F2(p a ,Pb)> with p a = p(J;l c ,l c + 51) and 
Pb = p{J + 5 J; l c , l c + 51), where for clarity we plotted only the 51 = 1,2,3 cases, as the other curves almost overlap 
the results for 51 — 3. Again, as in the previous case of the fidelity between the one-site states, here as well the sudden 
drop of the fidelity from its common value F2(p a , Pb) — 1 can be observed precisely in the vicinity of the critical point 
Jo = 1.9. We see it as a signature of both the change in the impurity one-site partial states, as well as the change of 
the total correlations between the impurity and the remaining sites. From our study of the one-site fidelity, we have 
seen that the bulk partial states are all approximately the same, irrespectively of the site and the parameter value, 
while only the impurity and first neighbor one-site states exhibit a visible dramatic change around the point of the 
quantum phase transition Jo. 

On the other hand, for every value of J we can consider the quantity 

r(TT + STi i +m = W(J; lc + 5l),p(J + 8 J; l c , lc + 51)) 

+ oj,t c , i c + oi) _ Fi( .. (J; ^ KJ + §J _ lc))Fi{p{J . lc + 5l)j KJ + SJ . lc + Sl)) > U»i 
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FIG. 3: One-site fidelity, given for fixed values of the parameter J, between any site and: a site in the bulk (upper figures); a 
site in the impurity (lower figures). Left figures are taken for J~ = 1.5 < Jo, figures on the right for J + = 2.5 > Jo, where 
Jo = 1.9 is the point of the quantum phase transition. Note the drop of the fidelity between the impurity and the bulk state, 
from F{~ ~ 0.96 for J~ — 1.5 < Jo, to f\ + ~ 0.65 for J + = 2.5 > Jo — another signature of the quantum phase transition. 



which can be seen as an indicator of total correlations, as in the case of uncorrelated two-site composite systems the 
above quotient reduces identically to 1. Therefore, the two-site composite system is necessarily correlated whenever 
F2 y 1. Note that C2 is not the measure of correlations — the two-site composite system can still exhibit nonzero 
correlations between its one-site subsystems even for F2 = 1. This is clearly seen in Fig. where C2 = 1 for all 
values of J except for J = Jo, although is was shown in 3J] that the considered two-site systems are correlated 
for every value of J. Yet, precisely around the point Jo of the phase transition, the quantity C2 exhibits an abrupt 
change, thus showing an enhanced change in the amount of the total correlations between the one-site subsystems of 
an overall two-site system, see Fig. [5] We also see that its numerical value around J = Jq approaches 1 as the step 51 
increases — ie the change of total correlations decreases as the distance between the two sites increases. This behavior 
was already observed before in the study of total and quantum correlations (entanglement) in this model [34| . In 
connection to this, we note here that the fidelity induces various distance measures on the set of quantum states 
(for an overview, see for example [22| and references therein), which can further induce fidelity-based multi-particle 
entanglement measures, based on the quantum state distinguishability, similarly to the case of the well known relative 
entropy of entanglement measures (36j . 

In the case of two-site subsystems with both one-site electrons located in the bulk, the fidelity taken between states 
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FIG. 4: Two-site fidelity between two states taken for two close parameter points J and J + SJ (SJ — 0.05), when one site is at 
the impurity and other in the neighborhood of it. Note the sudden drop of the fidelity around the point of the quantum phase 
transition, for the case of impurity + neighboring site. 



obtained for two close parameter values J and J + 6 J, is essentially equal to one (with SJ = 0.05). Analogously to the 
case of the fidelity between one-site partial states, here as well we see that the impurity does not affect considerably the 
partial states of the bulk of the lattice, even around the point of the quantum phase transition. Yet, we have observed 
a small deviation around Jo from the otherwise common value Fi ~ 1 which indicates that for larger subsystems 
(when considering the fidelities F3, F±, etc. of the 3-electron and 4-electron subsystems respectively), the impurity 
might start to increasingly influence the collective behavior of the bulk. 

Finally, as in the previous case, we have evaluated the Uhlmann mixed state geometric phase <f) g , again obtaining 
the trivial value <p g ~ everywhere, apart for the case when one electron site is located at the impurity: the structure 
of the two-site partial states' eigenvectors is not considerably affected by the quantum phase transition. 



CONCLUSION 



In this paper we have applied the concept of fidelity to approach quantum phase transitions, taking as the physical 
system a conventional 2-dimensional superconductor with a classical magnetic impurity placed at its centre. Unlike 
previous studies, where the overall states of the system were considered, here we have studied the fidelity for partial 
one- and two-site states. We have shown that the point of the quantum phase transition is still clearly marked by 
the sudden drop of the mixed state fidelity, even in the simplest case of the one-site partial states, as long as they 
are close to the impurity itself. The same holds for the two-site partial states — the fidelity exhibits a sudden drop 
in the vicinity of the quantum phase transition. Here, for the transition to be inferred it is necessary that one of the 
subsystem electrons lies at the impurity, or its neighborhood. This is expected to occur for this model since, as shown 
earlier, the impurity acts mainly as a local perturbation (but see also the behavior of the non vanishing two-site "long 
distance" entanglement measures discussed in [HI). We have also evaluated the Uhlmann mixed state geometric 
phase for the cases of both one- and two-site partial states and found it to be trivial (ie equal to one) for every value 
of the quantum phase transition driving parameter J, except when an electron site is at the impurity, thus showing 
that the structure of the subsystem eigenvectors has little relevance to the magnetic transition discussed. Finally, we 
have shown that the fidelity approach can also be useful in studying the total correlations between subsystems of the 
whole physical system, and how the point of quantum phase transition can be inferred from the change in the amount 
of total correlations (given by the quantity C2 (fT8j) ). which is closely related to the fidelity between the partial states 
discussed in this work. 
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FIG. 5: The C2 function between two states taken for two close parameter points J and J + 8 J (SJ = 0.05), when one site is 
at the impurity and other in the neighborhood of it. Note its sudden drop around the point of the quantum phase transition, 
as well as the decrease of its deviation from 1 as the other site gets further from the impurity. 



We believe that this work offers yet another confirmation of the usefulness of the fidelity approach to study quantum 
phase transitions. It extends this approach to partial states, showing that, in some cases at least, quantum phase 
transitions can be signaled out by the behavior of the mixed state fidelity taken between local partial states. It is an 
interesting and challenging task to try to find a unified view to this new approach and possibly connect it to some 
previously well studied thermodynamic quantities. Also, the generalization to thermal phase transitions is another 
possible path of further research. Finally, we see our simplified study of the total correlations in a composite quantum 
state, based on the fidelity, as just a first step towards exploring the possibility of using the fidelity-based quantum 
state distinguishability for quantifying correlations, both classical as well as quantum (entanglement), a subject of 
research that has on its own found rich applications within the field of condensed matter and many-body physics. 
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